Ground state energy of unitary fermion gas with the Thomson Problem approach 
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The dimensionless universal coefficient £ defines the ratio of the unitary fermions energy density to 
that for the ideal non-interacting ones in the non-relativistic limit with T = 0. The classical Thom- 
son Problem is taken as a nonperturbative quantum many-body arm to address the ground state 
energy including the low energy nonlinear quantum fluctuation/correlation effects. With the rela- 
tivistic Dirac continuum field theory formalism, the concise expression for the energy density func- 
tional of the strongly interacting limit fermions at both finite temperature and density is obtained. 
Analytically, the universal factor is calculated to be £ = g. The energy gap is A — y^fc/ /(2m). 
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With the further developments of the Bardeen-Cooper- 
SchriefTer theory, the possibility about the existence of 
the fermions superfluidity in the dilute gas system moti- 
vates widely theoretical studies and experimental efforts. 
Since DeMarco and Jin achieved the Fermi degeneracy [1], 
the ultra-cold fermion atoms gas has stirred intense inter- 
est about the fundamental Fermi-Dirac statistical physics 
in the strongly interacting limit. 

Across the Feshbach resonance regime, the interaction 
changes from weakly to strongly attractive according to 
the magnitude of the magnetic field. At the midpoint 
of this crossover unitary regime from Bardeen-Cooper- 
Schrieffcr(BCS) to Bose-Einstein condensation(BEC), 
the scattering length will diverge due to the existence of 
a zero-energy bound state for the two-body system. In 
this limit, the only dimensionful parameter is the Fermi 
momentum kf at T — 0. The corresponding energy scale 
is the Fermi kinetic energy £/ = k 2 /(2m), while m is the 
fermion mass. According to dimensional analysis, the 
system details do not contribute to the thermodynamics 
properties, i.e., the thermodynamics properties are uni- 
versal [2-19]. The energy density should be proportional 
to that of a free Fermi gas E/V = £ (E/V) free = £§ne/. 

This fundamental dimensionless universal con- 
stant £ has attracted much attention theoreti- 
cally/experimentally in recent years. Various theoretical 
approaches have been tried and the results differ from 
each other remarkably, for example, see Ref. [2-19] and 
references therein. The theoretical calculations are about 
£ ~ 0.3 — 0.6. Intriguingly, the experimental results 
are also quite different from each other, for example, 
£ « 0.74 ± 0.07[20], £ re 0.7[21], £ = 0.32±g;^[22], 
£ = 0.51±0.04[23], £ = 0.46 ±0.05 [24], £ = 0.46±£?|[25]. 
The recent lattice result is £ = 0.25 ± 0.03[26]. 

How to approach the exact value of £ analytically is 
a bewitching topic in the Fermi-Dirac statistical physics. 
To attack this intriguing topic is a seriously difficult prob- 
lem in many-body theory. The essential task is how to in- 
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corporate the nonlinear quantum fluctuation/correlation 
effects into the thermodynamics by going beyond any 
naive loop diagram expansions or the lowest order mean 
field theory. To our knowledge, the hitherto considera- 
tions looking for £ have been solely limited in the non- 
relativistic frameworks and with quite different results. 
How about a relativistic Dirac phenomenology attempt? 

Motivation : Essentially, the unitary physics with in- 
finite scattering lengths is quite similar to the universal 
strongly instantaneous Coulomb correlation thermody- 
namics in a compact nuclear confinement environment 
resulting from the competition of long and short range 
forces [27]. At the crossover point, the cross-section be- 
tween the two-body particle is limited by a ~ At: /k 2 [k is 
the relative wavevector of the colliding particles), while 
the gauge vector boson propagator is A M „ ~ g^v/k 2 in 
the gauge field theory. The former is short range but 
exact long range infrared correlation while the latter is 
long range one. The analogism motivates us to use the 
latter to model the former. 

In this context, it is also instructive to recall the in- 
termediate vector boson (IVB) hypothesis in weak inter- 
action theory. The local intermediate vector boson the- 
ory is related with the current-current (CC) contact in- 
teraction version through the corresponding connection: 
9w/ m w = I n the low energy limit k — ► 0, the two 
IVB and CC theories are identical. To model the unitary 
limit, we "let" g^y/ra^ —> oo with an arbitrary large 
charge gw or with an arbitrary small mass gap mw- 

It is usually assumed that there exists a uniform op- 
posite charged background to ensure the stability in 
discussing the electrons system thermodynamics. This 
is the well-known classical Thomson Problem[28]. Of 
course, there is not a finished standard answer for it- 
self/generalized version for over a century [29, 30]. Fur- 
thermore, a consummate method to gauge the infrared 
singularity in the dense and hot gauge theory is still to 
be looked for. However, we will turn the unitary limit 
topic into the same involved infrared problem. Then, the 
Thomson Problem is used as a potential quantum many- 
body nonperturbative arm to attack the infrared one. 
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To achieve the intriguing physics with the unusual in- 
mcdium relativistic Lorentz invariance breaking at uni- 
tarity, an infrared correlating "QED" Lagrangian is pro- 
posed in this Letter. Let the fermion have an "electric" 
charge g in addition to other internal global U(l) symme- 
try quantum numbers. According to the general stabil- 
ity principle, the system should be stabilized by a Active 
uniform opposite charged background in the meantime. 
This particular assumption makes it possible for us to 
deal with the challenging infinity. 

The natural units c = h = ks = 1 are used. 

To perform the path integral presentation as a nonper- 
turbative starting point, the considered effective actions 
involve the interaction of Dirac fermions with an auxil- 
iary Proca-like Lorentz vector boson field [31, 32] 



^matter = ^(h^ ~ TTl)^; 

c Af 1 
l 



4 ^ ' 



r — — 2 
J-i,a-^ - 2 m background^ 



A U A» + A^, (1) 



where m is the bare fermion mass. The A^ is the vector 
field with the field stress 



d„A v - d u A 



(2) 



In the action Ci t A-i>, the electric vector current con- 
tributed by the fermions is 



(3) 



Based on the local gauge invariant free Lagrangian, the 
many-body interactions can be introduced with a hidden 
local symmetry (HLS) manner[33, 34] 

Li = -\f^F^ + \D^H\ 2 + V(H) + A„ J". (4) 

The many-body Lorentz violation environment mod- 
ulates the bare two-body interaction. The vector boson 
mode mass gap ^background * s n0 ^ an m t erna l degree of 
freedom. It appears as a free parameter and controls the 
interaction strength between the particles. It indicates 
the interaction of the opposite electric charged back- 
ground or the stochastic many-body potential caused by 
the strongly fluctuating/correlating effects. Below, the 
suffix "background" will be emitted for brevity just with 
a tilde symbol trig. The un-physical coupling constant 
g also represents the unitary limit/infrared characteristic 
"two-body" bare potential. 

The electric current conservation/gauge invariance is 
guaranteed by the Lorentz transversality condition with 
HLS formalism 



dpA" = 0, 



(5) 



which can be naturally realized by taking the relativistic 
Hartree instantaneous approximation (RHA). 



In terms of the functional path integral[27, 35-37], the 
auxiliary effective potential reads 



Sl/V 



-m 2 A 2 



2T [ \ln{l + e-^ Ek -^)) 

J k 



+ ln(l + e-^ fc +^) 



(6) 



where "2" represents the (hyperfinc-)spin projection of 
the fermions with J k = J d 3 k/(2ir) 3 and (3 = 1/T being 
the inverse temperature. The tadpole diagram with the 
boson self-energy for the full fermion propagator leads to 



a 9 
A = T n, 



(7) 



with the fermions (electric charge number) density de- 
fined by the thermodynamics relation dfl/d/i\A = — n 



2 / [/(//*, T) — /(//*, T)] . 

Jk 



(8) 



In the above expressions, 
1 



/(/AT) 



e P(E k -v*) _|_ I 



1 



e /3(E fc + M *) + 1' 



are the distribution functions for (anti-)particles with 
From Eq.(7), the effective chemical 



E k = Vk 2 + m 2 
potential /i* is defined with a gauge invariant manner 

„2 



H* = fj, + /ii = \x y-n, 



(9) 



where [i is the global chemical potential. The spirit 
is quite similar to the Kohn-Sham density functional 
theory [38]. 

Using Eq.(6) and with the thermodynamics relation 



1 d{m) 

one obtains the energy density functional 



(10) 



e = ~—n z + 
2m| 



2 f E k [f(jS,T) + f(jS,Tj\. (11) 
Jk 



The second term in Eq.(ll) appears as very much 
the analytical formalisms for the free Fcrmi-Dirac gas. 
However, the correlating effects are implicitly included 
through the effective chemical potential esp. for T^O. 

The mass gap parameter m 2 ^ is a Lagrange multi- 
plier that enforces relevant constraints and reflects quan- 
tum fluctuating effects consequently. The remaining cen- 
tral task is how to determine the unknown many-body 
stochastic potential characterized by the coupling con- 
stants m 2 ~ and g 2 . The answer can be found from the 
auxiliary effective potential, i.e., the composite system 
should be "charge neutralized" through the Active back- 
ground with the artificial conditional extremum [39] 



= 0. 



(12) 
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With Eq.(12), one can have 

J ' ' (2k 2 + m 2 ) 



m| 



9 



hi d|k| 



o 



Ek 



[/(/i*,T) + /(/i*,T)] 
(13) 



which is the negative of the gauge invariant Debye 
(Thomas- Fermi) screening mass squared m 2 D , i.e., there 
is a minus sign between ra|, and m 2 ^. The Debye mass 
parameter can be also directly calculated with the vector 
boson polarization through the Dyson-Schwinger equa- 
tion (relativistic random phase approximation- RPA) [40] 



with the Fermi kinetic energy £/ = fc 2 /(2m) and particle 
number density n — kj-/ (3-7T 2 ). Therefore, one can obtain 
the ground state binding energy (energy per particle) 



eb = to = 

n 



15 42 to 3 \m) 



£/ .(19) 



By keeping up to the lowest order oikf/m, the ratio of 
the binding energy to the Fermi kinetic energy is ^ . Fur- 
thermore, one can obtain the the universal dimcnsionless 
coefficient £, i.c, the ratio of the energy density to that 
of a free Fermi gas 



n^(p ,p) 
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Tr 



1 



r 



to 



m 



with the full fermion propagator Eq.(9) by noting 



m 



D 



n°°(o,| P |-o). 



(14) 



(15) 



In Eq.(14), the 0-component of the four- momentum k — 
(fcojk) in the fermion loop is related to temperature T 
and effective chemical potential [i* via fco = (2n+l)nTi+ 
fi* . It is very interesting that the Thomson stability con- 
dition can give the gauge invariant Debye mass, which is 
exactly consistent with the Dyson-Schwinger equation. 



At T — 0, the parameter m ~ is 

J Jd 



(16) 



k) 



with kf being the Fermi momentum and Ef - 

the relativistic Fermi kinetic energy. 

The Eq.(ll) with Eq.(13) (m| = -m 2 D ) arc our 
main results, from which one can further study the 
strongly correlating fermions thermodynamics in the uni- 
tary limit. The collective interaction contribution is neg- 
ative for the physical energy density functional. Espe- 
cially, there is not any remained parameter because the 
un-physical coupling constant g appears simultaneously 
in the denominator and numerator within the relevant 
analytical expressions through m 2 ^. From Eq.(16), one 
can see that the magnitude of this fraction ratio m 2 ^/g 2 
characterize the density of states and consequently the 
fluctuating contributions. 

With the mathematically well defined energy density 
functional Eq.(ll) and with Eq.(13), we now return to 
the final discussions. At T — > 0, the energy density is 



k 5 



(2k 2 +m 2 )k f E f - m Hn k -^ 



18ir 2 E f 



8tt 2 



.(17) 



In the non-relativistic limit kf/m <C 1, one can expand 
e according to the Taylor series of kf/m 



—/( — 



free 



4 5 _ 4 
15 X 3 ~ 9' 



(20) 



This analytical result is in the range of the existed 
theoretical ones £ ~ 0.3 — 0.6[4-19]. It is exactly 
consistent with that of the quantum Monte Carlo 
calculations [5] and in reasonable agreement with the 
updating experiments [22-25]. We also note a simi- 
lar £d=oc = | was obtained by Steele within the ef- 
fective theory long ago, with D being the space-time 
dimensions [41]. 

At unitarity, the energy gap should be also propor- 
tional to the Fermi kinetic energy according to the di- 
mensional analysis. It can be empiristicly derived from 
the total energy density with the usual odd-even stag- 
gering in the thermodynamics limit[5, 11] A = yg£/, 
with e f being the Fermi kinetic energy. The critical tem- 
perature is T c i=a 0.157TV approximated with the BCS 
universal relation T c = e 7 A/7r, where 7 is the Euler 
constant. The T c is in agreement with the updating 
theoretical[19, 42]/experimental[23, 43] results. It should 
be pointed out that the differences for the energy gap 
A / critical BCS phase transition temperature T c can be as 
large as several times in the literature. In the meantime, 
the validity of the weak coupling BCS relation between 
A and T c deserves to be further studied theoretically. 

Let us further analyze the pressure by keeping up only 
to the lowest order of kf 



P = 



1 



30m 



(21) 



The power index of n is also consistent with the ex- 
isted theoretical calculations [44]. One can see that the 
strongly correlating effect as well as the fermion mass 
affects the coefficient a significantly. 

The sound speed squared for the ideal Fermi gas is 



via 



1 dP-, 



FG 



dn 



1 2 

ft' 



(22) 



with the Fermi velocity vf — kf/m. In the unitary limit, 
the pressure is P = ^Pfg, from which one can find the 
sound speed is reduced remarkably 



4 5 4 1 4 

e = ran + — net + —n— - -n^ 
15 1 42 to 3 m 



(18) 



VFG 



V2 



v f > 0. 



(23) 
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It has been argued that the sound speed squared might 
be negative due to the theoretical spinodal instability in 
the unitary limit [4]. Although the sound speed is signif- 
icantly reduced due to the correlating effects, it is still a 
real number which indicates the system stability. 

The interior correlations have been refreshingly incor- 
porated in the thermodynamics as an external source 
manner indirectly through a mirror Thomson back- 
ground. In other words, we have derived the effective 
interaction strength, with which we have given the gen- 
eral but concise analytical thermodynamic expressions 
at both finite temperature and density. This analytic 
method can be easily extended to the near unitary limit 
regime [45]. 



In conclusion, there is a similarity between the strongly 
interacting ultra-cold atoms physics in the scattering 
length limit \a\ — > oo and the infinite instantaneous 
Coulomb interaction in a compact confinement environ- 
ment due to the competition of long range forces with 
short range ones. The classical Thomson Problem as a 
potential quantum many-body arm redounds to address- 
ing the universal thermodynamics, with which the key 
coefficient £ = | and the energy gap A = yg£/ are ob- 
tained in the QED-likc framework with HLS formalism. 
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